Outbreaks of Hantavirus induced by seasonality 
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Abstract 

Using a model for rodent population dynamics, we study outbreaks of Hantavirus infection 
induced by the alternation of seasons. Neither season by itself satisfies the environmental re- 
quirements for propagation of the disease. This result can be explained in terms of the seasonal 
interruption of the relaxation process of the mouse population toward equilibrium, and may shed 
light on the reported connection between climate variations and outbreaks of the disease. 
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I. INTRODUCTION 



Hantaviruses are rodent-borne zoonotic agents that may cause diseases in humans such 
as hemorrhagic fever with renal syndrome and hantavirus pulmonary syndrome 
Hantaviruses have been identified at an increasing rate in recent years, and as of now about 
thirty different ones have been discovered throughout the world. One of these, the Sin 
Nombre virus, was not isolated until 1993 after an outbreak in the Four Corners Region in 
the USA DD . The host of this particularly dangerous virus is the deer mouse, Peromyscus 
maniculatus, the most numerous mammal in North America. The virus produces a chronic 
infection in the mouse population, but it is not lethal to them. It is believed that transmission 
in the rodent population is horizontal and due to fights, and that the subsequent infection 
of humans, where the mortality rate can be as high as fifty percent, is produced by their 
contact with the excreta of infected mice. Moreover, so far there is no vaccine or effective 
drug to prevent or treat the hantavirus pulmonary syndrome. Therefore, a major effort has 
been made to understand the population dynamics of deer mice colonies in order to design 
effective prevention policies jlj]. 

It has been noted that environmental conditions arc directly connected to outbreaks of 
Hanta Q]. For instance, the 1993 and 1998 outbreaks occurring in the Four Corners Region 
have been associated with the so-called El Nino Southern Oscillation Q. Related to this 
and other such observations, the effects of seasonality in ecological systems have been a 

-year oscillations of mammal populations prev- 
predator seasonal dynamics [7], and persistence of parasites in plants between seasons |8J] 
are examples that illustrate the importance of seasonality in population dynamics. 

Recently, Abramson and Kenkre have proposed a phenomenological model for mice pop- 

rn 

ulation that successfully reproduces some features of Hantavirus propagation [9( . In partic- 
ular, that model explores the relation between resources in the medium, carrying capacity, 
and the spread of the infection in the rodent colony. Herein we study the effects of seasonal- 
ity in that model. Our motivation is not only to provide more realism to the model, but also 
to investigate the counterintuitive effects that dynamic alternation may cause in a biologi- 
cal system. Brownian motors jl(| and switching-induced morphogenesis Q| are examples 
that show that alternation in time of "uninteresting" dynamics may produce "interesting" 
outcomes. Along these lines, we will show that alternation of seasons, neither of which by 
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itself fulfills the environmental requirements on the carrying capacity for spreading of the 
infection, may produce an outbreak of the disease. The mechanism driving this behavior 
is the interruption of the relaxation process that equilibrates the mouse population from 
season to season: if the duration of seasons becomes shorter than the relaxation time of the 
population, the disease spreads. 

The paper is organized as follows. In Sec. |H] we briefly review the model for mouse 
populations introduced in Je| . In Sec. IIHI we explain how seasonality is introduced in that 
model and analyze the conditions for Hanta outbreaks to take place due to the alternation 
of seasons. The exact solution of the model and a particular example that illustrates the 
phenomenology is given in Sec. II VI Finally, in Sec|VJ we summarize the main conclusions 
and propose some directions for future work. 



II. THE MODEL 
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The model introduced in |9| for the temporal evolution of a population of mice subjected 
to the Hantavirus infection reads: 

^ = 6 M-cM s -^-„M s M„ (la) 
at K 

dMi „, MjM 

__L = - cMl - -JL- + aMsMj, (lb) 

where Ms and Mj are the population densities of susceptible and infected mice respectively, 
M = Ms + Mi is the total population of mice, and a, b, c and K are positive constants. 
The terms on the right hand sides of Eqs. (fTajl and (jib)) take into account the following 
processes: births with rate constant b, depletion by death with rate constant c, competition 
for the resources in the medium characterized by the carrying capacity K, and transmission 
of the infection with rate coefficient a. It is worth noting that infected pregnant mice 
produce Hanta antibodies that keep their fetus free from the infection; that is, all mice are 
born susceptible 1], as indicated by the absence of a birth term in Eq. (jlb|) . Note also the 
absence of a recovery term in the model since, as mentioned earlier, mice become chronically 
infected by the virus. 

The system of equations (jla)) and (jlbj) has four equilibrium points. Two of them are 
irrelevant for the analysis: the null state Mj = Ms = 0, which is always unstable if b > c (a 
condition that we will assume throughout this paper), and a meaningless state with Mi < 
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FIG. 1: Stable equilibrium population of susceptible (solid line) and infected (dashed line) mice 
as a function of the resources present in the medium, K. The values of the parameters are the 
same as those used in Ref. a=0.1, 6=1, and c=0.5. The value of the critical carrying capacity 
is K c = 20. 

for any value of the parameters. The other two equilibria are 

M s = K(b - c) Mi = 0, (2) 

M s = h - Mj = K(b-c)-^ } (3) 

The stability of the equilibrium points (j2J and (JHJ) depends on the value of the carrying 
capacity. If K < K c = - (rz~) then (J2J) is stable and © unstable. If K > K c it is the other 
way around. That is, when the available resources, K, are below the critical value, K c , 
the infection does not propagate in the colony, the whole population of mice grows healthy, 
and its size increases proportionally with those resources. As soon as K surpasses K c the 
virus spreads in the colony, the susceptible mouse population saturates, and the fraction of 
infected mice becomes larger as K increases (see Fig. 
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III. SEASONAL ALTERNATION 



The Four Corners Region, where an important number of cases of Hantavirus Pulmonary 
Syndrome have occurred, has a desert climate. The largest climate variations within this 
region come from the alternation between dry and rainy seasons. We will therefore assume 
alternation in time of these two seasons. It is important to remark that a two-season 
assumption is not crucial, and that the analysis with four seasons is also straightforward 
within the formalism introduced herein. During each of the two seasons under consideration 
we assume there to be no climate variations, so that each season can be characterized by a 
set of time-independent parameters {p^ = {oj, 6j,Cj, where i — 1,2. We implement 

square-periodic season alternation where the duration of each season is T/2. Again, this 
particular alternation pattern is not essential for the mechanism. Other schemes of season 
alternation, e.g. different duration of the seasons or random switching between seasons, do 
not qualitatively change the phenomenology. 

Any quantity p(t) alternating in the way described above can be written as: 

p(t) = p + + p-n{t), (4) 
where p± = \{pi ± p 2 ) and p{t) is a periodic square wave 



p{t) 



1 : < t < f 
-1 : \ < t < T 



Let us now suppose the following conditions for the sets {pi} according to seasonality: 

ai < a 2 , h < b 2 , ci < c 2 , K 1 > K 2l (5) 

where 1 stands for the rainy season and 2 for the dry one. The biological motivation for 
these conditions is the following. The dry season provides less resources for the colony than 
the rainy season [K 2 < Ki), and as a consequence the death rate is higher (c 2 > Ci), and the 
transmission rate is also larger since fights for the available resources are expected to increase 
(a 2 > ai). However, notice that we consider the birth rate to be larger during the dry season 
(b 2 > bi). There are two reasons for this. First is the assumption that the colony makes 
an attempt to maintain its population. Second is an implementation of the maturation 
process in the model: baby mice do not contribute to the propagation of the disease jlj], 
and, therefore, even if births were actually more numerous during the rainy season, the 
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FIG. 2: (Ki, i^-region of disease propagation due to seasonal alternation. The values of the other 
relevant parameters are given in Eqs. © and (JJJ). Every value of K\ or K<i within the enclosed 
area by itself does not satisfy the required environmental conditions that support spreading of the 
virus in the mice colony. Yet, those same points lead to outbreaks of Hanta when alternation of 
seasons is sufficiently fast. 

contribution of this fact to the propagation of the infection will only be important in the 
next (dry) season, when mice have matured and are ready to fight. It will be shown later 
that these assumptions lead to a situation with high population during the rainy season and 
low population during the dry one, in agreement with the available data 

Moreover, we assume that K\ < mm(K cl ,K c2 ), i.e., the resources are all times (during 
both seasons) below the minimum critical threshold that triggers the propagation of the 
disease. We will show that nevertheless it is possible for the infection to spread. 

The equilibrium populations of the susceptible and infected mice are determined by the 
set of values {pi}- When switching from season to season, the populations evolve trying 
to reach a new equilibrium. Therefore, the dynamics is driven by the competition between 
two characteristic times. On the one hand there is an external time scale determined by the 
seasonal forcing, t e = T/2. On the other hand, the relaxation toward equilibrium after a 
switching of seasons involves a relaxation time. The latter measures the time required for 
the mouse colony to relax to the equilibrium state associated with {pj} after having been 
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driven during the previous season by the conditions {pi}, that is, t r (i — > j) where i,j = 1, 2. 
The internal time scale is defined as the fastest relaxation process, i.e., ti = min[£ r (z — > j)]. 
"Fast" or "slow" seasonal alternation then refers to the comparison between these two time 
scales. If t e 3> ti the mouse population has enough time to accommodate to the new 
conditions from season to season and relax to equilibrium. Moreover, since we have imposed 
the condition that the resources at any time of the year are below the critical thresholds 
K C i, there will be no infected mice. In the other limit, t e <C ti, seasonal changes occur too 
fast, the relaxation process is interrupted, and no equilibrium can be reached from season 
to season. Then note that an adiabatic elimination can be implemented 3|, and fj,(t) in 
Eq. (@J) can be replaced by its average value, = 0. Therefore, in the limit of fast season 

alternation the system is driven by the set of averaged values {p+} = {a + , b + , c + , (1/ ' K) + }, 
and the critical carrying capacity is given by K c+ = a+ (b+-c+) • ^ s a consec Luence, h is 
possible to find regions of parameters where K c+ is smaller than the effective value of the 
carrying capacity associated with the averaged values: 



i -i 



1 1 
2\K 1 + K~ 2 



2K1K2 



and the infection propagates. 

General conditions leading to this behavior can be posed, but the expressions are rather 
cumbersome. We prefer, for the sake of simplicity, to show a particular typical case. We use 
the following values for the parameters: 

oi = ^ h = !, ci = g, (6) 
73 

a 2 = 4, b 2 = — , c 2 = 1, (7) 

that lead to K c \ = 6 and K c2 = 73/4 respectively. The dynamics are completely determined 
once the value of the carrying capacity during each season is specified. According to the 
previous discussion, these parameters can be chosen such that the following conditions hold: 

[(1/K) + ] 1 > K c+ , K x > K 2 , K x < mm(K cl , K c2 ). 

These conditions lead to the points (K±, K 2 ) that fulfill the seasonal requirements given by 
Eq. so that slow alternation of seasons leads to infection- free states while fast alternation 
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leads to Hanta outbreaks. This region is plotted in Fig. El Notice in particular that the 
point {K\ — 4, K% — 1) lies inside the region and that Ki < K C i. In the next section we 
illustrate the seasonality-induced propagation of the disease for this particular point. 



IV. THE CRITICAL PERIOD 

So far we have determined that outbreaks of Hanta induced entirely by seasonal changes 
can occur if the duration of the seasons are short enough. Now we establish the meaning 
of "short enough" quantitatively. Since K c+ is strictly smaller than the effective value of 
the carrying capacity, there should be a finite value of T c such that for any T < T c the 
population of infected mice is greater than zero, but for periods above this critical period 
the infected population goes to zero. 

In order to obtain the value of the critical period we need to solve the system of equations 
()la|) and (jlbj) . In spite of its nonlinearities the system can be solved analytically by means 
of a reciprocal transformation ^| and the following exact solution is obtained: 

M ' (i ^ s ° iW) = s^rShwr (8a) 

M s (t, /„, S„; M) = SM ° ei [ , (8b) 

where Jo and Sq are the initial conditions for Mi and Ms respectively, and the following 
definitions have been introduced, 

aK-l 



Q(t) = e- ct (M (e& - l) + S 
S = K(b-c), M = I + S . 



Because the external forcing due to the alternation of seasons is periodic, Floquet theory 
ensures the existence of a periodic solution. The values of Iq and Sq compatible with the non- 
equilibrium periodic solution can be obtained by evolving the system during the first half of 
a period under dynamics 1 and the second half under dynamics 2, and forcing periodity on 
the solutions after a whole period of evolution, that is, 

Mi (~ Mi I , S ; { Pl }) , M s J , Sq- { Pl }^j ; {p 2 ^j = I , (9a) 
M s Mi C^, Jo, S ; { Pl }\ , M s C^, J , S ; { Pl }) ; {p 2 }\ = S . (9b) 



FIG. 3: Population of susceptible (top) and infected (bottom) mice versus the period of the 
seasons. The dashed and continuous lines indicate the populations at the end of the rainy and dry 
seasons respectively. The critical period for which the virus begins to spread due to seasonality is 
T c ~ 12. The values of the relevant parameters are a\ = 1/4, b\ = l,c± = 1/3, K\ = 4 and a 2 = 4, 
&2 = 73/72, C2 = 1, K 2 = 1 for the rainy and dry seasons respectively. 

In order to close the system in the non-equilibrium steady state Mi (t, T; {pi^}) and 
M s (t,T; {pi^}), the values of I and So that solve that system of equations (J9~aj) and (J9b"j) 
must be then re-introduced in Eqs. ()8a|) and (j8b|) . The critical period is then the largest 
value of T satisfying the condition 

Mi(t,T;{p 1>2 }) >0. 

We illustrate the procedure with the example mentioned above where the parameters are 
given by Eqs. © and (J7J), and with K\= A. and K 2 = 1. The results are shown in Figs. H3 
and |3J The values of Mi and Ms as a function of the period of the seasons are depicted in 
Fig. El where the populations of the susceptible and infected mice at the end of each season 
are given. As seen in that figure, the value of the critical period is T c ~ 12. Note that if the 
alternation is slow, T > T c , all mice grow healthy. On the other hand, if the alternation is 
faster than the relaxation time required by the colony to accommodate its population from 
season to season, T < T c , the virus spreads and Mi > 0. Notice that in the limit T —>■ the 



FIG. 4: Population of susceptible (top) and infected (bottom) mice versus time for a period of 
evolution. The values of the relevant parameters are the same as those used in Fig. |21 The critical 
period is T c ~ 12. Left panel: Results for a very long period (T = 40). Right panel: Very short 
period (T = 0.1). Central panel : Results for a near-critical period (T = 10). 

dynamics is driven by and the populations of susceptible and infected mice are given 
by Eq. (JHJ) with a = a + , b = b + , c = c+, and 1/K = {1/K) + . Let us stress again that the 
carrying capacity is below its critical threshold at any time. 

In Fig. |Uwe plot, for different period lengths, the exact solutions Mj (t,T; {pi^}) and 
Ms (t, T; {pi^}) as a function of time through one period of evolution. The first semi-period 
corresponds to the rainy season and the second to the dry season. When seasons last long 
(left panel), there are no infected mice and the susceptible population simply oscillates 
between the two equilibrium states given by Eq. (J2J). For sufficiently short seasons (right 
panel), there is propagation of the disease and the values of Mj and Ms fluctuate around 
the equilibrium points determined by Eq. (jSJ) and the set of parameters {p+}. Finally, when 
the period of the seasons is near, but below, the critical period (central panel), the infected 
population is small and the population Mj oscillates in a more pronounced manner. 

V. CONCLUSIONS 

By introducing seasonality in a paradigmatic model for Hantavirus propagation in mice 
colonies, we have shown that the alternation of seasons may cause outbreaks of the disease. 
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The striking feature of that behavior lays in the fact that neither season satisfies the con- 
ditions for the infection to spread in terms of the availability of resources. The mechanism 
responsible for the phenomenon is the competition between two time scales: an external 
one, the duration of a season, and an internal one, the relaxation time for the mouse colony 
to equilibrate its population from season to season. We have shown that if the duration of 
the seasons is longer than the relaxation time, no propagation of Hantavirus occurs. On the 
other hand, if the relaxation process is interrupted by a fast seasonal alternation, the dis- 
ease spreads. We have analyzed the general conditions for which the phenomenon occurrs. 
Moreover, we have illustrated the mechanism with a particular example that can be solved 
exactly. 

This work may help to clarify the reported relation between climate and propagation 
of Hanta in mice populations. However, to elucidate whether the proposed phenomenon 
actually takes place in nature we depend on data that unfortunately are not available in 
the literature. One can envision further modifications of the model that may improve its 
features, such as, for example, the inclusion of spatial dependence or of noisy contributions 
to the dynamics. Finally, we stress that the general idea underlying the mechanism is model- 
insensitive and can therefore be extended to other systems where seasonality plays a relevant 
role. Work along these directions is in progress. 
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